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Abstract 

A density disturbance in a system of gravitating mass, induced by a moving se¬ 
lected body gives rise to a dynamical screening of Newtonian potential of this body. 
When applied to the solar planetary system it means that as a result of the motion of 
the Sun in the Galaxy its effective force potential appears more weak than the New¬ 
tonian potential. The relevant modifications of main relations of the solar dynamics 
are considered here and it is found in particular that the reestimated period of the 
Earth revolution around the Sun rises in 1 second per year and semimajor axis of the 
Earth orbit increases on 4 kilometers. Similar relations are obtained for other planets 
too. It may be supposed that the inclusion of these effects can help to explain the 
observable anomalous acceleration of spacecrafts Pioneer 10 and 11. 


1 Introduction 

It is a matter of general experience that a resting test charge in a plasma produces the 
Debye screened potential. When this charge is moving the static screening decreases [1] 
and the potential becomes anisotropic [2]. These properties are due to collective effects 
in systems of charge particles. Here the similar effects for systems of gravitating particles 
are discussed. For simplicity we shall restrict our consideration to the model steady state 
system of a homogeneous Maxwellian gas of gravitating bodies. 

The problem of motion of the test star in a fluctuating force held of another stars was 
hrst studied by Chandrasekhar [3] and discussed later by Marochnik [4] using methods of 
kinetic gas theory as an initial problem for a distribution function describing the test star. 
In contrast to his approach, here we consider a perturbation of the equilibrium state of a 
gas of gravitating bodies as a response to a force held of the test moving body following to 
our general approach developed in refs. [5-7]. 

The main advantage of our approach is a derivation in Sec.2 of the expression for a renor¬ 
malized (ehective) potential of the Sun and calculation of its diherence from the Newtonian 
solar potential for planets of the solar system. 

Next we calculate in Sec.3 corresponding corrections to the periods and semimajor axis 
of revolutions of the Earth and other planets. 


2 Gravitational screening 

In contrast to plasma or electrolyte where the screening of the Coulomb potential of a 
charge particle is presumably due to the presence of particles of an opposite charge and 
takes a place both for the moving test particle (dynamical screening) and for the hxed one 
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(static Debye-Hiickel screening) the screening of the Newtonian potential of a selected test 
body in a system of gravitating bodies is possible for the moving test body only. 

In the model discussed below the Galaxy is considered (neglecting the dark matter) as 
an equilibrium gravitating system of stars with a chaotic distribution of velocities described 
by the Maxwellian /o(v) = po(27rh^)“^/^ exp{—, where po is the density of stellar 
matter and v velocities dispersion. The Sun is excluded from the whole system and 
considered as the additional gravitating body of mass Mq inserted at the instant to at 
the point tq = 0 with the determined steady velocity u . Introducing the test particle 
gives rise to a perturbation /i(r, v,t) of the system, so its distribution function takes the 
form /(r, v,t) =/o(v) +/i(r, v,t) . 

Then, in the linear approximation in the perturbation the system is described by the 
Vlasov equation 


<9/i(r,v,t) (9/i(r,v,t) (9/o(r,v,t) _ 

dt dr dv 

This is supplemented by the Poisson equation for the effective potential 

<h = —AnGm J d^vfi{r, v, t) — 4:7rGMQ6{r — u (t — to))- 


( 1 ) 

( 2 ) 


We are interested in a steady-state behavior of the gravitating system. It means physically 
that processes are considered at times t > r , where r = 1/V47rGpo is the characteristic 
relaxation time. For our Galaxy (where po = 1.9 ■ 10“^^gcm“^ ) it is of order 10^ years 
what falls far short of the time of existence of the solar system (in the opposite case 
transition processes would be to take into consideration). To exclude transition processes, 
we put to —oo . Then the formal solution to equations (1) and (2) is 

^ = 4:7rGm f (f V f dt''V^{r — ^{t — t'),t') ■— 4:7rGMQ6{r — ut) (3) 

J J —OO UV 


or in the Fourier transform: 


[k‘^ — inGm 


dh 


k ■ dfo/dY 
u — kv + iu 


]$(k,o;) = STT^GM^diu 


ku). 


(4) 


whence 


Here 


<h(k, u) 


Stt^GMq 

k^eg{k,u!) 




ku). 


eg(k,uj) 


1 47rG /• 3 k■^fo/^^r 

1 “ -TT / ^ ^ - 1 — 

k^ J oj — kv + lu 


(5) 

( 6 ) 


is the gravitational permittivity of the system. This function determines a reply of the 
gravitating system on a gravitational disturbance. The concept of the gravitational per¬ 
mittivity is well-known in the theory of gravitating system (see, e.g. [4, 8, 9]). Its zeroes 
dehne a dispersion equation of a gravitating system of which complex roots determine a 
spectrum of elementary excitations. 

An inhnitesimal imaginary value iu is added to the frequency u in equations (4) 
and (6), which may be considered as a result of introducing a ’’shadow” of the Boltzmann 
collision integral —z//i into the r.h.s. of the Vlasov equation. Such addition assures a 
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selection of retarded solutions and accounts for Landau damping in collisionless plasma. 
Then the permittivity ^^(k, cu) becomes a complex function and can be expressed as 

where kj = {AnGpo/v'^)^^'^ and W{z) = 1 + i{7TY^‘^zexp{—z^} erk{—iz) . 

Calculating integrals of reverse Fourier transform 

$(r,t) = r du [ ~ $(k,o;) (8) 

J—oo J 

we choose the axis 2 : along u and introduce the dimentionless variables Z = {z — 
ut)K,j, X = XKj, K = k/^j, V = u/{\/2v) . Then, taking into account cylindrical sym¬ 
metry we get 


(7) 


$(r,t) = 


GMqKj f 

27r2 J- 

GMqKj 
(X2 + Z2)5 L 


d(j) J dO j dKK^ sin 6 


exp{iK (X sin 9 cos (j) + Z cos 9 )} 


X2 - W{Vcos9) 

l + (A'= + Z")3/(A',Z,l/)l, (9) 


where 


I{X,Z,V) = ^ r d<i) r d9 f dKW{V cos 9) sin 9 
J-n Jo J 


exp{iK (X sin 9 cos (j) + Z cos 9 )} 
X2 - W{Vcos9) ■ 


( 10 ) 

This integral determines a deviation of the effective potential from the Newtonian one 
described by the hrst term in the square bracket of the equation (9). After taking the 
Cauchy integral over K we hnd 


I{X,Z,V) 


^ r d(t) r d9 sin 9V^ 

47r2 J-TT Jo 


-2i sinh(Av/^)Ci(-iAy^) + 


+cosh(Av/^) (7r + 2iShi(Av/^))], 


( 11 ) 


where A = X sin 9 cos 4> + Z cos 9 . To calculate numerically values of this integral over 
an arbitrary range of variations of the parameters V and Z should present no problems 
but for the discussed subject of planetary motion in the solar system we can conhne our 
consideration to concrete values of these parameters. 

For stellar (barionic) matter in the Galaxy we have kj ~ 2.6 ■ 10“^^ cm“^ , thus the 
values Z and X can be estimated as 10“® < jZj < 10“® , 10“® < |X| < 10“® for planets 
of the solar system of which orbits remote from the Sun at distances 10^^ < r < 10^® cm . 
Since a peculiar velocity of the Sun in the Galaxy is 1.95 ■ 10®cm/c , and the velocity 
dispersion in the Galaxy is h ~ 1.55 ■ 10®cm/c , we put V = 1 . 

As a result of calculation of the integral (11) at X = 1 we get / = —0.3677 . When 
varying V within the range from 0.3 to 3.0, the numerical value of the integral I varies 
from -0.3 to -0.47. 


3 Correction to revolution periods 

According to the equation (9) and obtained numerical estimation of the integral I the 
energy of the renormalized gravitational interaction between the Earth and the Sun can be 
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represented as 


U{R) =Uo + SU, Uo = -A/r, 5U = A 7 , A = GMqMq,, 7 = -Jkj. (12) 


Thus, in such approximation the correction 6U to the Newtonian energy of interaction is 
constant value independent on r . 

Before proceeding to a calculation of its influence on the Kepler period of the Earth 
revolution it is necessary to check whether it disturb the closure of the hnite trajectory of 
the Earth around the Sun. 

In the classical two-body problem (after its reduction to the problem of the one body 
motion in a central force held) an angle of precession of the perihelion Aip (see, e.g. [ 10 ]) 
is dehned as 


A(p = 2 


4dr 


^2m{E-Uo-SU)-^' 


(13) 


where 

Mq+Mq 


^ max S.ncl T min 

reduced mass. 


are the maximum and minimum of the radius-vector r , 


m = 


M = mr'^ip, E 



2mr‘^ 


+ U. 


The last two values, the moment M and energy E , are the motion integrals. 

The condition of closure of a trajectory is that the angle of precision to be of the form 
Acp = 27rk/n where k,n are integers. For Newtonian potential this condition is fulhlled 
( A(p = 271). 

As a result of expanding of the integrand of the equation (13) in SU the zeroth member 
gives 27r , and the hrst member determines an additional precession of the perihelion as 


5(p = 


d /2m 


dM 




(14) 


where the expression within the brackets (...)o is estimated for the unperturbed motion. 
With the use of the parameter p and eccentricity e of an elliptic closed orbit we have 


Then 


r = 


P 

(1 + e cos 99 ) ’ 


P 


m 2 _ r 2EM2 
mA ’ ^ V mA'^ 


6(p 


221A (r __) 

mAdM \ Jo {l + ecoscpyj 

27 d M^tt 2771 d ( 
mA 9M^(l-e2)3/2^ ^ 


0 . 


(15) 


Therefore, the perturbed trajectory remains closed. 

We estimate now a change of the revolution period T under influence of the pertur¬ 
bation to the energy SU . It is not difficult to check that all the mathematical treatment 
leading [10] to the third Kepler law for undisturbed Newtonian potential 


Tn = ttA, 


m 


2\En 


En — 


mr 


M 2 

2mr‘^ 


+ Uo{r) 


( 16 ) 
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remains valid when we change Eq for E = Eq + 6U and leads to the modified form 


T 


71 A, 


m 


2\Eo + 6U\^ 



SU\ 


- 3/2 




To 



Ml] 

2lEolJ' 


( 17 ) 


whence 


It can be noted that ao 
we get 


T - To _ 3yl7 
To 2|To| 


(18) 


A 

2 |So| 


is the semimajor axis. Substituting here the value of Oq 


T-Tq 

To 


= 3700 = 4.25-10" 


■ / ■ 

r G 1 

1/2 

r Po 1 

[-0.367?] 

[6.67- lO-sJ 


o 

1 

to 

CO 


1/2 


X 


X 


V 

-1 

ao 

[l.55 ■ 106J 


[l.496 ■ 1013J 


(19) 


which corresponds to AT^ ~ 1.3 sec/year (1 year ~ 3.15 ■ 10^ sec) 

Thus, accounts for screening of the Newtonian potential gives rise to lengthen of the ET 
year on 1 second approximately what compensate (in the order of the value) the systematic 
difference AT = UT — ET = UT — TAI ~ 0.8 sec. 

Another consequence of the renormalization of the Newtonian potential is modihcation 
of parameters of the Earth orbit. In particular, the Earth semimajor axis a® oc |T|“^ 
increase on the value 

Aa® Ay 


iTn 


~ 2.83 ■ 10" 


( 20 ) 


that is, Aa® ~ 4.25 ■ 10^ cm. 

Similar calculations may be performed for other planets of the solar system. Their 
results are represented in the Table. 


Planet 

AT/T 

AT, sec 

Aa/a 

Aa, km 

Venus 

3.1 ■ 10-3 

0.6 

2.1 - 10-3 

2.2 

Earth 

4.3 ■ 10-3 

1.3 

2.8 - 10-3 

4.2 

Mars 

6.5-10-3 

3.8 

4.3 - 10-3 

9.9 

Jupiter 

2.2 ■ 10-^ 

0.83-10^ 

1.5 - 10-1 

1.1 - 10^ 

Saturn 

4.1 ■ 10-^ 

3.8 ■ 10^ 

2.7-10-1 

3.9-10^ 

Uranus 

8.2 ■ 10-1 

2.2 ■ 103 

5.4-10-1 

1.6 - 103 

Neptune 

1.3 ■ 10-6 

6.7-103 

8.5-10-1 

3.8-103 


Thus, we see that relative time and space corrections to the far giant planets motion 
resulted from the renormalization of the solar potential are of order 10“® . It may be 
supposed that they can be detected with the use of spacecrafts on far solar orbits. 

The above results are based on some approximations. The most important of them are 
our choice of galactic disk volume density po = 1.9 ■ 10“^^ g cm“^ having an influence on the 
value of the Jeans wave number kj and total ignoring of a dark matter. As it was shown in 
[6], inclusion of the dark matter gives rise to a sufficient decrease of an effective Jeans wave 
number and enhancing of the dynamical screening. Therefore, an experimental checking 
of the above results could give an additional information related to the dark matter in the 
Galaxy. 
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It would be interesting to apply our results to discussion of the observed [11, 12] anoma¬ 
lous Pioneer 10 and 11 acceleration with a magnitude ~ 8 x 10“® cms“^ . Its relative value 
is of order of 10“^ , that is much greater our corrections to orbits of the giant planets 
being at the same distance ( 20 30 AU) from the Sun. An expected value of the space¬ 

craft acceleration was calculated with the use of two independent programs, that are JPL’s 
Orbit Determination Program (ODP) and Aerospace Corporation’s Compact High Accu¬ 
racy Satellite Motion Program (CHASMP). Both programs use the planetary efemeris, and 
timing inputs. Without knowledge of the programs we can only suppose that taking into 
account of our corrections to these data can produce a change of the expected acceleration 
and, probably, provide a reasonable explanation for the observed anomalous acceleration. 
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